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Abstract 

The averaging problem in cosmology is of considerable importance for the correct 
interpretation of cosmological data. A rigorous mathematical definition of averaging in 
a cosmological model is necessary. In general, a spacetime is completely characterized 
by its scalar curvature invariants, and this suggests a particular spacetime averaging 
scheme based entirely on scalars. We clearly identify the problems of averaging in a 
cosmological model. We then present a precise definition of a cosmological model, and 
based upon this definition, we propose an averaging scheme in terms of scalar curvature 
invariants. This scheme is illustrated in a simple static spherically symmetric perfect 
fluid cosmological spacetime, where the averaging scales are clearly identified. 
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1 The averaging problem in cosmology 



Cosmological observations [H [2], based on the assumption of a spatially homogeneous and 
isotropic Friedmann-Lemaitre-Robertson- Walker (FLRW) model plus small perturbations, 
are usually interpreted as implying that there exists dark energy, the spatial geometry 
is flat, and that there is currently an accelerated expansion giving rise to the so-called 
ACDM-concordance model. Although the concordance model is quite remarkable, it does 
not convincingly fit all data (3J. Unfortunately, if the underlying cosmological model is not 
a perturbation of an exact flat FLRW solution, the conventional data analysis and their 
interpretation is not necessarily valid. For example, the standard analysis of type la super- 
novae (SNIa) and CMB data in FLRW models cannot be applied directly when backreaction 
effects are present, because of the different dynamical behaviour of the spatial curvature 
[I]. Indeed, supernovae data can be explained without dark energy in inhomogeneous mod- 
els, where the full effects of general relativity (GR) come into play. For example, it has 
been indeed shown that the Lemaitre-Tolman-Bondi (LTB) solution can be used to fit the 
observed data without the need of dark energy [5], although it is necessary to place the ob- 
server at the center of a rather large-scale underdensity. Therefore, the averaging problem in 
cosmology is of considerable importance for the correct interpretation of cosmological data. 
The correct governing equations on cosmological scales are obtained by averaging the Ein- 
stein field equations (EFE) of GR (plus a theory of photon propagation; i.e., information on 
what trajectories actual particles follow). By assuming spatial homogeneity and isotropy on 
the largest scales, the inhomogeneities affect the dynamics though correction (backreaction) 
terms, which can lead to behaviour qualitatively and quantitatively different from the FLRW 
models. 
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1.1 General Approaches 



The gravitational field equations on large scales are obtained by averaging the EFE of GR. It 
is necessary to use an exact covariant approach which gives a prescription for the correlation 
functions that emerge in an averaging of the full tensorial EFE. The Universe is not isotropic 
or spatially homogeneous on local scales. An averaging of inhomogeneous spacetimes on 
large scales can lead to important effects. For example, on cosmological scales the dynamical 
behavior can differ from that in the spatially homogeneous and isotropic FLRW model [6]; in 
particular, the expansion rate may be significantly affected. Consequently, a solution of the 
averaging problem is of considerable importance for the correct interpretation of cosmological 
data. The averaging problem in GR and cosmology is of fundamental importance. 

There are a number of approaches to the averaging problem [HI E] • In the approach of 
Buchert [8] a 3+1 cosmological space-time splitting is employed and only scalar quantities 
are averaged. The perturbative approach involves averaging the perturbed Einstein equa- 
tions; however, a perturbation analysis cannot provide any information about an averaged 
geometry. In the space-time or space volume averaging approach tensors, and in some cases 
only scalar quantities, are averaged; this procedure is not generally covariant hence the re- 
sults are somewhat limited and the conclusions unreliable. In all of these approaches, in 
analogy with Lorentz's approach to electrodynamics, an averaging of the Einstein equations 
is performed to obtain the averaged field equations. 

However, to date the macroscopic gravity (MG) approach is the only approach to the 
averaging problem in GR [7] which gives a prescription for the correlation functions which 
emerge in an averaging of the non-linear field equations (without which the averaging of 
the Einstein equations simply amount to definitions of the new averaged terms). The MG 
approach is a fully covariant, gauge independent and exact method. The space-time av- 
eraging procedure adopted in MG for any differentiable manifold is based on the concept 
of Lie-dragging of averaging regions, and it has been proven to exist on an arbitrary Rie- 
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mannian space-times with well-defined local averaged properties. Averaging of the structure 
equations for the geometry of GR brings about the structure equations for the averaged 
(macroscopic) geometry and the definitions and the properties of the correlation tensors. 
The averaged Einstein equations can always be written in the form of the Einstein equations 
for the macroscopic metric tensor when the correlation terms are moved to the right-hand 
side of the averaged Einstein equations to serve as the geometric modification to the averaged 
(macroscopic) matter energy-momentum tensor [TJ. 

The formal mathematical issues of averaging tensors on a differential manifold has re- 
cently been revisted However, we note that integrating scalars on spacetime regions is 
always well-defined. 

1.2 Scales 

In any theory of physics, the scales over which the physical theory is applicable must be 
specified [TO]. There is a hierarchy of cosmological scales of physical interest. Consequently, 
we must specify the cosmological scale over which averaging occurs (i.e., we must specify 
the averaging scale £ or averaging region, which then determines the type of averaging or 
smoothing that occurs). In particular, a given geometry may be inhomogeneous on both 
small scales £ s (i.e., local scales such as local density inhomogeneities) and large scales £[ 
with £ s « £i (where £\ ~ £jj, and £h is the Hubble scale). The averaging scale £ will satisfy 
£ s « £, but it is possible that £ « £i, so that the geometry is still inhomogeneous (but 
'smooth') on large scales. In cosmology it is assumed that the averaging scale £ is bigger 
than scale of largest observed structures (clusters of galaxies) and voids and that £ < £h 
[TT] . There is also the homogeneity scale, £h om , the largest scale any inhomogeneities are 
observed. It is usually assumed that £h om < £ (< £h)- 

The physical description of a cosmological model depends on the averaging scale [12J. 
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The scales £ s , £i, £ in the cosmological model and the range of validity (i.e., £ s « £ < £\) 
must be specified. Of course, the range of scales relevant to cosmology are the largest 
scales of averaging, larger than the largest scale of cosmological structures and comparable 
to a fraction of £u- In the context of perturbation theory, particular attention is paid to 
both the scale of inhomogeneities in the background and the scale of inhomogeneities of the 
perturbations [T3] . 

1.3 Spatial curvature 

In P3] the MG equations were explicitly solved in a FLRW background geometry and it 
was found that the correlation tensor (backreaction) is of the form of a spatial curvature. 
Thus, the averaged Einstein equations for a flat spatially homogeneous, isotropic macroscopic 
space-time geometry has the form of the Einstein equations of GR for a non-flat spatially 
homogeneous, isotropic space-time geometry. 

The relevance of spatial curvature in realistic models of the universe that describe the 
dynamics of structure formation since the epoch of last scattering was discussed in [13]; in 
particular, in arguments about spatial curvature in perturbation theory, the quasi- Newtonian 
approximation must be used with care since spatial curvature is an inherently relativistic 
phenomenon (that does not occur in Newtonian physics). We note that a spatially dependent 
constant spatial curvature k can alleviate the tension in observational data [15]. Indeed, if 
the spatial curvature parameter k is allowed to be a function of position, then considerable 
spatial curvature (locally) is permissable (consistent with CMB observations) [TB] . 
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1.4 Discussion 

Clearly, backreaction (averaging) effects are real, but their relative importance must be 
determined [17]. Observational data suggests a normalized spatial curvature |fifc| ~ 0.01 — 
0.02 (i.e., of about a percent). Combining these observations with large scale structure 
observations then puts stringent limits on the curvature parameter in the context of adiabatic 
ACDM models; however, these data analyses are very model- and prior-dependent [4J, and 
care is needed in the proper interpretation of the data. There is a heuristic argument that 
fi fe ~ 10~ 3 — 10~ 2 [IB1EI5], which is consistent with CMB observations [TJ[2J[2J|I] and agrees 
with estimates for intrinsic curvature fluctuations using realistically modelled clusters and 
voids in a Swiss-cheese model. It must be appreciated that such a value for fife, at the 1% 
level, is relatively large and may have a significant dynamical effect on the evolution of the 
universe and the interpretation of cosmological observations. 

Note that in a scenario in which ~ 0.01 — 0.02, the current contribution from the 
spatial curvature is much greater than the energy density of radiation, and is comparable 
to the energy density in luminous matter. In addition, such a value cannot be naturally 
explained by inflation. From standard analysis, depending on the initial conditions and the 
details of a specific model of inflation, |fi — 1| would be extremely small. Indeed, any value 
for fife at the 1 % level would be very difficult to explain within the theory of inflation; 
therefore, any non-zero residual curvature at this level can only be naturally explained in 
terms of an averaging effect. 

1.5 Null geodesies 

Ultimately we wish to determine the effects of averaging or backreaction on the evolution 
of the universe and the interpretation of cosmological observations. All deductions about 
cosmology are based on light paths. Only the redshift and the energy flux of light arriving 
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from a distant source are observed, rather than the expansion rate or the matter density. It 
is often assumed that intervening inhomogeneities average out. However, inhomogeneities 
affect curved null geodesies [181 E] and can drastically alter observed distances when they 
are a sizable fraction of the curvature radius. In the real universe, voids occupy a much 
larger region as compared to structures, hence light preferentially travels much more through 
underdense regions and the effects of inhomogeneities on luminosity distance are likely to be 
significant. 

The effect of averaging null geodesies in inhomogeneous models was discussed in [T7] . 
GR is treated as a microscopic (classical) theory. Real photons travel on null geodesies 
in the microscopic geometry. However, because all observations are of finite resolution, 
observations necessarily involve averages of measured quantities. Therefore, in interpreting 
real observations, it is necessary to model properties of (not only a single photon but of) a 
'narrow' beam or bundle of photons (i.e., a local congruence of null geodesies). From the 
geometric optics approximation we can obtain the optical scalar (Dyer-Roeder) equations 
that govern the propagation of the local shearing and expansion (of the cross-sectional area 
of the beam) with respect to the affine parameter along the congruence due to Ricci focussing 
and Weyl tidal focussing |19J . Since the nonlinear optical scalar equations require integration 
along the beam, the optics for a lumpy distribution does not average and there may be 
important resulting effects. Similar issues have been discussed recently from a different point 
of view p^8l [20] . The formulation of the EFE on the null cone was discussed in [21]. Clearly 
averaging can have an important effect on photon propagation and hence observations. 

In particular, assuming GR to be a microscopic theory on small scales with local metric 
field g (the micro-geometry) and matter fields, a photon follows a null geodesic k in the 
local geometry, and after averaging we obtain [17] a smoothed out macroscopic geometry 
(with macroscopic metric (g)) and macroscopic matter fields, valid on larger scales. But, in 
general, the "averaged" vector (k) need not be null, need not be geodesic (and even if it is, 
need not be aflinely parametrized) in the macro-geometry. 
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2 Scalar curvature invariants and averaging 



In [22], it was shown that the class of four-dimensional (4D) Lorentzian manifolds that cannot 
be completely characterized by the scalar polynomial curvature invariants constructed from 
the Riemann tensor and its covariant derivatives must be of Kundt form. This implies that, 
in general, a spacetime is completely characterized by its scalar curvature invariants and 
this suggests a particular spacetime averaging scheme based entirely on scalars. Let us first 
review the main mathematical background. 

For any given spacetime (M, g) we define the set of all scalar invariants 

T — J T) E> PA"' n rniva.fi p p/xi/»/3;7 p p/^a/3;7<5 1 (n i \ 

_£. =\n } n^ v n , O^a^O , Itf U , a ^--yIl , lXfiuafi-^S 1 ^ ■>■■■}■ 

Consider a spacetime (M, g) with a set of invariants X. Then, if there does not exist a 
continuous metric deformation of g having the same set of invariants as g, we will call the 
set of invariants non- degenerate. Furthermore, the spacetime metric, g, will be called X- 
non- degenerate. This implies that for a metric which is X- non-degenerate the invariants 
characterize the spacetime uniquely, at least locally, in the space of (Lorentzian) metrics. 
This means that these metrics are characterized by their curvature invariants and therefore 
we can distinguish such metrics using their invariants. 

It was proven, on a case-by-case (depending on the algebraic type; characterized by their 
Petrov and Segre types or, equivalently, in terms of their Ricci, Weyl (and Riemann) types 
|23j). that a 4D spacetime is either X-non-degenerate, which implies that the spacetime 
metric is determined locally by its curvature invariants, or the metric is a Kundt metric. 
This is a striking result because it tells us that metrics not determined by their curvature 
invariants must be of Kundt form. These Kundt metrics therefore correspond to degenerate 
metrics in the sense that many such spacetimes can have identical invariants. The Kundt 
class is defined by those metrics admitting a null vector £ that is geodesic, expansion-free, 
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shear-free and twist-free 

Let us also consider the 'inverse' question: given a set of scalar polynomial invariants, 
what can we say about the underlying spacetime? In practice, it is somewhat tedious and a 
lengthy ordeal to determine the spacetime from the set of invariants. In 4D we can partially 
characterize the Petrov type in terms of scalar curvature invariants. In most circumstances 
we only need some partial results or necessary conditions or we are dealing with special cases. 
For example, we found that if 27 J 2 7^ J 3 , or if 27 J 2 = J 3 but the differential invariants S\ 7^ 
or S2 7^ 0, then the spacetime is X-non-degenerate. Having determined when a spacetime is 
completely characterized by its scalar curvature invariants, it is also of interest to determine 
the minimal set of such invariants needed for this classification. It is also of interest to study 
when a spacetime can be explicitly constructed from scalar curvature invariants. 

Let us return to the question of averaging. We have noted that a general spacetime is 
completely characterized by its scalar curvature invariants. Since we know how to average 
scalar quantities, we can average all of the scalar curvature invariants that can then represent 
an averaged spacetime (with that set of averaged scalar invariants). We shall return to this 
in Section 4. However, we note that cosmological models belong to the set of spacetimes 
completely characterized by their scalar curvature invariants, suggesting that we can average 
a cosmological model using scalar invariants. We must first define a cosmological model 
rigorously, which we do next. 



1 We recall that in the Riemannian case a manifold is always locally characterized by its scalar polynomial 
invariants. 

2 We note that this exceptional property of the the degenerate Kundt metrics essentially follows from the 
fact that they do not define a unique timelike curvature operator. 
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3 Cosmological models and averaging 



There are a number of technical problems with averaging of tensor fields on a differentiab 
manifold M. Clearly we need a covariant approach. Is the average of the metric the metric? 
There is the question of averaging verses smoothing. We want to avoid issues with respect to 
coordinates; i.e., if we choose more smooth coordinates, we are not averaging but just using 
different coordinates to represent the same geometry. It may be possible to avoid several 
of these technical problems by adopting an approach based on scalar curvature invariants. 
Note that although such an approach may work for any differentiable manifold which is 
1 -non- degenerate, we shall focus on the cosmological problem for the most part here. 

Therefore, we wish to discuss the averaging problem in the context of cosmology. Al- 
though this context may simplify the issue of averaging in technical terms, there are some 
new problems of principle that are then introduced. First, the cosmological model is a 
mixed model, in that the matter is already assumed to be averaged, but the geometry is 
not (necessarily) |j. Therefore, we need a consistent model for the matter, represented on 
the characteristic averaging scale, and its appropriate (averaged) physical properties []. It is 
known that the separation between the gravitational field and the matter is not scale invari- 
ant and the notion of a perfect fluid is not scale invariant [25]; averaging (in the presence 
of a gravitational field) modifies the equation of state of the matter [SB]. 111 addition, since 
averaging does not conserve geodesies (in fact, averaging need not even conserve the metric 
signature), we need further assumptions in order to be able to compare the models with 
observational data. 



3 Note that if the average of the metric is not the metric itself, then the conformal structure of the macro- 
and micro- geometries are not the same and photons will not follow null geodesies in the macrogeometry in 
general. 

4 Since the time dependence is smooth, perhaps some smoothing of the microgeometry has already been 
assumed. 

5 Note that a discrete (non-continuous fluid) model for the cosmological matter has been discussed recently 

EH. 
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First, a precise definition of a cosmological model is necessary; i.e., a framework in which 
to do averaging. This definition includes an appropriate way to do averaging and how to 
deal with photons. In particular, it is necessary to make all of the assumptions in the model 
clear. 

The precise definition of a cosmological model we shall adopt, based in part on [TBI \T7\ I2T] . 
is given by the following conditions Cl — C5 (cf. [10J): 

3.1 Definition of a cosmological model 
Cl. Spacetime Geometry: 

The spacetime geometry (M, g) is defined by a smooth Lorentzian metric g (characterizing 
the macroscopic gravitational field) defined on a smooth differentiable manifold M. 

The macroscopic metric geometry is obtained by an appropriate spacetime averaging of 
the microgeometry; thus part of the definition of a cosmological model consists of specifying 
the averaging scheme (which must be consistent with the physical assumptions of the model 
encapsulated in the conditions C3 and C4 below) and the cosmological scale over which 
averaging or the smoothing occurs (i.e., we must specify the averaging scale £ or averaging 
region). In particular, a given microgeometry may be inhomogeneous on both small scales £ s 
and large scales t\. The averaging scale £ will satisfy £ s « £, but it is possible that £ « £i, 
so that the macrogeometry is still inhomogeneous (but 'smooth') on large scales. As noted 
earlier, the scales for which the cosmological model is defined must be explicitly specified. 
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C2. Timelike Congruence: 

There exists a timelike congruence (u) (in principle locally, but by definition we can extend 
this to the whole manifold), representing a family of fundamental observers. 

This congruence is associated with the 4-velocity of the averaged matter in the model; 
i.e., the matter admits a formulation in terms of an averaged matter content which defines 
an average (macroscopic) timelike congruence. If there is more than one matter component 
giving rise to more than one macroscopic timelike congruence we identify a fundamental 
macroscopic timelike congruence. There is always one, which has a physical meaning. This 
leads to a covariant 3 + 1 split of spacetime [10]. Mathematically this means that the 
spacetime is X -non- degenerate and hence the spacetime is uniquely characterized by its 
scalar curvature invariants. 

In addition to the formal parts CI and C2 of the definition of a cosmological model 
(M, g, £, u), we must also specify the physical relationship (interaction) between the macro- 
scopic geometry and the matter fields, including how the matter responds to the macroscopic 
geometry. 

C3. Macroscopic field equations: 

There exists an appropriate set of macroscopic field equations relating the averaged matter 
and appropriately averaged (or macroscopic) geometry. This is based on an underlying 
microscopic theory of gravity (such as, for example, GR), and an appropriate formalism to 
average the geometry and find corrections (correlations) due to averaging to the Einstein 
tensor in the resulting field equations: 

G\ + C\ = T% (3.2) 
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in the usual geometrized units in which c = 1 = 8nG/c 2 , where G a b = R a b — \b~\R and R a b 
is the Ricci tensor of the averaged macrogeometry, C a b is the correlation tensor, and T a b is 
the energy momentum tensor (already assumed averaged). 

This is the emphasis of the analysis here. Note that in this context only the Ricci tensor 
needs to be averaged. 

Discussion 

The energy-momentum tensor that appears in the EFE (defined formally by the variation 
of the matter action with respect to the metric) depends on the gravitational field. Indeed, 
in GR the gravitational energy contributes to the total energy-momentum. So in this sense 



some aspects of the averaged gravitational field are already included in the model. 



In general, it does not follow from the contracted Bianchi identities that the energy- 
momentum is conserved with respect to the macrometric: T c b . c ^ 0. Is it possible to define 
a new effective energy-momentum tensor T a b which is conserved? For example, T a b = G a b = 
T a b — C\ (which satifies the macro-conservation equations with respect to the macrogeometry 
by definition). However, it must be ensured that the relationship between the effective matter 
and the (averaged) macrogeometry is consistent with the underlying microphysical model 
and the averaging scheme. At a heuristic level, this could be (within the cosmological model 
under consideration) modelled qualitatively through an appropriate (new type of) effective 
equation of state. 

In particular, the separation between the gravitational field and the matter is not scale 
invariant. Indeed, the notion of a perfect fluid is not scale invariant [25], and in the presence 
of a gravitational field the equation of state of the matter is consequnetly in general scale 
dependent (26]. Thus, averaging affects the effective equation of state. In the simplest 
cosmological models we can write T a b = diag[—p,p,p,p] and C a b = diag[— p c ,p c ,p c ,p c ], and 



To minimize this effect, the mixed components of T a b are usually considered. 
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we have that T a b = diag[p — p c ,p—p c ,p—p c ,p—p c }. A physical equation of state relates p and 
p, whereas the effective equation state is a relationship between p = p — p c and p = p — p c . [] 

In this heuristic framework all of the effects of averaging go into the redefined energy- 
momentum tensor T a b and the effective equation state of the macro-matter (subject to any 
changes in the equations of motion of the macro-matter, although we again note that T\ 
is conserved relative to the macrogeometery) . Therefore, in this reinterpretation, any non- 
standard equation of state is not due to any exotic matter, but is the effect of averaging 
(cf. [28]). In particular, an equation of state renormalization at the level of about 1% 
~ 0(£/£h) naturally fits into this new averaging interpretation (and especially an additional 
C a b corresponding to spatial curvature arising from averaging [23J). 



C4. Equations of motion: 



matter moves (and also 
In principle, averaging 



We also need to know the trajectories along which the cosmologica 
the light trajectories, which determine observational relations). |® 
the 4- velocities of the microscopic matter particles does not necessarily give the macroscopic 
4-velocity of the cosmological matter [] (and the average motion of a photon need not be a 
null geodesic in the averaged geometry). 

Discussion 



7 If we split up T into a standard matter part plus an 'additional' matter term, in general there exists an 

interaction between them (i.e., each matter term is not necessarily conserved |27j ) . 

8 This may add additional geometrical quantities (other than the curvature) that need to be considered 

(and need to be averaged). Perhaps this can be done using appropriate scalars, but this is the subject of 
future investigation. 

9 In GR a free point particle moves on a timelike or null geodesic. But a system may behave as point 
particles on small scales but not on larger scales (alternatively, the definition of a point particle is also not 
scale invariant). Therefore, macroscopic (averaged) matter need not move on geodesies of the macrogeometry. 
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This needs to be consistent with the averaged field equations: if the matter satisfies the 
conservation equations in the macroscopic geometry (which follows from the averaged field 
equations), then the equations of motion must be consistent with this (i.e., cosmological 
dust follows macroscopic timelike geodesies). Note that G a b satisfies the contracted Bianchi 
identities (with respect to the macrogeometry) . 

The fundamental congruence is essentially the average of the timelike congruences along 
which particles move in the microgeometry. This implies a relationship between matter and 
macrogeometry. If the particles move on timelike geodesies in the microgeomerty, does the 
fundamental congruence consist of timelike geodesies (of the macrogeometry). Does (a beam 
of) photons move on null geodesies of the macrogeometry. This raises the question of whether 
it is possible to perform an appropriate averaging so that on average the macroscopic matter 
moves on timelike/null geodesies of the macrogeometry. 

Ultimately we need a theory for the propagation of photons. The motion of light in the 
macroscopic geometry must be consistent with the limiting motion of timelike particles in the 
macrogeometry. Alternatively, we could average the Einstein-Maxwell equations and take 
the geometric optics limit [19]; i.e., we assume that the actual photons satisfy the Einstein- 
Maxwell equations in the optical limit, and thus follows a null geodesic in the macrogeometry. 
We must then average the Einstein-Maxwell equations [26] and obtain a suitable form for the 
macroscopic Einstein-Maxwell equations and an appropriate (corresponding) optical limit in 
the macroscopic regime. 

C5. Observations: 

Finally, we need to be able to relate averaged quantities with physical observables, which 
ultimately must be consistent with cosmological data. 

10 Since T is conserved the motion of matter is restricted; we still need a model for the motion of photons. 
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3.2 Covariance 



We have advocated an approach to averaging a cosmological model utilizing scalar invariants. 
Thus far we have focussed attention on averaging the EFE (e.g., the Ricci tensor). We shall 
briefly address the more general question of averaging or smoothing a differentiable manifold 
in the next section. However, the second aspect of averaging in a cosmological model involves 
the equations of motion of matter and the geodesic structure. 

Since there exists a global timelike congruence u, global covariance is broken (and the 
spacetime admits of a 1 + 3 split). This timelike congruence is defined physically, and in 
principle is measurable. The resulting theory need only be invariant with respect to re- 
stricted covariance |29j. Perhaps, for describing the equations of motion of particles and null 
rays, we only need to consider scalars invariant with respect to restricted sets of coordinate 
transformations (this would correspond to the kinematic variables such as the shear and 
expansion scalars associated with the macroscopic covariant derivative of u [10J). 

This suggests that a description of the geodesic structure might be possible utilizing 
averaged the restricted (kinematical) scalar quantities. We shall not pursue this further here 
(cf. 0). 

3.3 Example: FLRW models 

In the standard FLRW model there are a number of simplifications and assumptions. How- 
ever, the cosmological model is not fully defined in the sense above and hence has a lack of 
physical predictability. 

The past approaches to averaging have been ideally suited to the FLRW models (with 
small, vanishing in the limit, perturbations). Therefore, an extension of these approaches 
to cosmological models other than FLRW models is not possible [2S|. Also, some of these 
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approaches deal exclusively with dust. Therefore, in these approaches to averaging and 
cosmology it is impossible to obtain the potentially most important effects. 

Let us examine the assumptions CI - C5 in the usual FLRW plus perturbations model. 
The macrometric g is the FLRW metric (CI) and u also has a geometric meaning (C2). In 
the usual point of view there are no correlations due to averaging (i.e., the correlation tensor 
is zero; C\ = 0) or, more precisely, they are negligible (C3). In this case it follows from the 
contracted Bianchi identities that energy-momentum is conserved: T c b . c = 0, which relates 
the matter to the averaged geometry. All other effects are assumed negligible (C4). 

However, there is no formal proof that such assumptions arise from a rigorous averaging 
scheme of some appropriate (physically motivated) microgeometry. In addition, there are 
some important effects which are not necessarily small perturbations. In [23] it was argued 
that important effects at 1% level cannot be neglected. 

Since there is no scale included in the model, in a sense the model is incomplete. Indeed, 
the model does not even have the ability to determine whether there is a scale above which 
the geometry is exactly FLRW or whether at all scales the geometry is only approximately 
FLRW (with a given perturbation scale). Indeed, in cosmological perturbation theory, both 
the scale of the background and the scale of perturbations need to be specified, neither of 
which are given since no notion of order of approximation is included [TO f [T3" | [T2] . In addition, 
the motion of photons must be independently postulated. 

Regarding C5 we ask the question: does this model agree with observations? If it does 
not, then even if the model agrees in some approximate sense with most observations, there 
is no structure within which to discuss the potential small discrepancies with observed data, 
which is a deficiency of the model. If the model does, then it would be remarkable, although 
there is still the need for a physical explanation of the dark energy. 

In the standard approach, there might be a spatial curvature term on the left-hand side 



17 



of the macroscopic field equations with intrinsic curvature parameter k, and an effective 
spatial curvature term on the right-hand side of the macroscopic field equations (due to an 
effective equation of state arising from averaging) with curvature parameter k (which could 
be thought of as a 'renormalization' to the intrinsic curvature parameter). If observations 
indicate that fij k 1 - 2% [U [2J, then there is no physical mechanism within the model 
(particularly if there is an inflationary period) to produce an intrinsic curvature parameter 
k of this magnitude, whereas an effective curvature parameter k of about a percent arises 
naturally from averaging [T5] . 
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4 An approach to averaging cosmological models using 
scalar invariants 



Since for a cosmological model the geometry is completely characterized by its scalar cur- 
vature invariants [22], we shall adopt the approach that we shall only average the scalar 
curvature invariants (thereby avoiding the problems mentioned earlier). Therefore, we have 
a microgeometry completely characterized by its set of scalar curvature invariants X. We then 
average these microgeometry scalar curvature invariants to obtain a new set of macrogeom- 
etry scalar curvature invariants X, which now completely characterizes the macrogeometry. 
In general, this macrogeometry is unique (for a given averaging region). Let us consider 
first the more general mathematical question of averaging the geometry; we shall then focus 
attention on the physical case of a cosmological model. 



4.1 Averaging the geometry 

In the general mathematical context we want to describe the averaged geometry (represented 
by the Riemann tensor and its covariant derivatives) and interpret the results. Let us 
consider, X = {R, R^B? , . . .}, which is an ordered set of functions on M. X is then also 
an ordered set of functions. If we write X = {R, R^ V R V ' . . .}, then of course it does not 
follow that, for example, R^vR 11 ^ = R fl u R u ^ (under averaging). But the question is: does the 
ordered set of functions X correspond to the associated scalar curvature invariants for some 
metric g (which will then serve to define the macrometric g). This raises some interesting 
mathematical questions (see below). In principle, this is not necessarily true in general. 

Indeed, we note that it is not generally true that the set of averaged scalar invariants X 
actually determines a Riemannian geometry. In addition, if A = B, then (A) = (B). |HJ If 



For convenience we use the simpler notation (A) = A here and in the last section. 
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A 2 = B 2 , then (A 2 ) = (B 2 ), but this does not imply (A) 2 = (B) 2 . In this particular case it 
is true because if A, B are scalars, then A 2 = B 2 implies A = ±B and (A) = ± (B); but if 
(as in the static example below) (R2) 2 = (-R1) 3 , while it is true that ((R2) 2 } = ((Rl) 3 }, it 
does not follow in general that (R2) 2 = (-R1) 3 . However, relationships of this form will be 
true for many of the applications of interest (see the next section). 

Indeed, it is plausible that the ordered set of functions X correspond to the associated 
scalar curvature invariants for some macrometric g for the class of X-non-degenerate geome- 
tries that constitute the class of cosmological models defined here. Since the geometries 
are X-non-degenerate and in 4D properties of the geometry can be represented in terms of 
scalars, and since relations between different terms (functions) in the set X (e.g., R and R 2 
are functionally dependent) and the corresponding terms in the set X (e.g., R and R 2 ) are 
functionally related in exactly the same way and syzygies (e.g., describing the algebraic type) 
are maintained under averaging, it follows that in general the set X uniquely gives rise to a 
macrometric g (which will have similar algebraic properties to the micrometric g). 

We note that a subset of 4D geometries are uniquely determined by their curvature only 
(i.e., no covariant derivatives are necessary), so that we only need to consider the scalar 
polynomial invariants constructed from the curvature tensor (i.e., an appropriate subset of 
X) o For algebraically general geometries this is indeed true, and it is certainly true for 
many of the important applications we are interested in. In addition, as noted above, we are 
primarily interested in the Ricci curvature (and the 4 Ricci scalar invariants) in applications 
in cosmology. 



12 This has not been proven here, but it is plausible; if there are any exotic counterexamples, then the 
definition of a cosmological model, and the rules of averaging, could be suitably amended to avoid this. It 
is true for the cosmological models of physical interest. An alternative approach is suggested in the scalar 

averaging procedure described below. 

13 That is, we do not need to consider invariants like diWeyl = C a b c d;eC abcd '' e or diRicci = R a i re R ah]e , 

which depend on derivatives of the metric greater than 2, which may add new problems regarding averaging. 
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If a spacetime is not X-non-degenerate it can not be uniquely determined by its scalar 
polynomial invariants. There is one subclass of X-non-degenerate spacetimes which, in a 
certain sense, can be determined by its scalar polynomial invariants; namely the type D k 
spacetimes (i.e., spacetimes for which the curvature and all of its covariant derivatives are si- 
multaneously of algebraic type D) [22]; this subclass includes the static spherically symmetric 
spacetimes of particular interest here. 



4.1.1 Uniqueness 

The macrogeometry is completely specified by the set X. Suppose that X (or some appro- 
priate subset of X) is derivable from a macrometric g. Although this macrogeometry is 
unique, it is possible for two different (non-isomorphic) microgeometries to give rise to the 
same macrogeometry. However, the algebraic properties (i.e., the algebraic properties of 
the Riemann tensor and its covariant derivatives) of the underlying microgeometry and the 
resulting macrogeometry must be the same. 

Since in the cosmological models under consideration the microgeometry is completely 
characterized by the set of scalar curvature invariants X, the algebraic properties of the un- 
derlying microgeometry is generally determined through the syzygies of X. o Any syzygy 
S is maintained under averaging (integration), giving rise to a corresponding syzygy S of 
the set X and then the averaged geometry will obey a similar set of (averaged versions) 
of these syzygies; that is, the averaged geometry will be of the same type as the microge- 
ometry, and hence the macrogeometry is at least as algebraically special as its underlying 
microgeometry] 15 ! In addition, if certain terms in X satisfy certain algebraic conditions due 
to the differential Bianchi identities, then the corresponding terms in the set X would sat- 



14 A particular cosmological spacetime (an exact solution or a subset of spacetimes with arbitrary functions 

in 4D) is often completely determined by syzygies involving scalar invariants. 

15 In principle, new syzygies may appear via averaging and so it is possible, in principle, for the macroge- 
ometry to be more algebraically special than its underlying microgeometry. 
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isfy a corresponding set of algebraic conditions, consistent with the Bianchi identities of the 
macrogeometry. Therefore, in applications in cosmology this approach will produce a unique 
macrogeometry. 



4.1.2 Proposal: Scalar Averaging Procedure 



Let us consider the ordered set of functions X on M in the form: 



X = {-R, i?2, i?3, R 2 , R iL v R v ^ . . . , C 2 , . . .}. 

First, let us omit any scalars from this set that are not algebraically independent (e.g., 
{R 2 , R^R? ' . . .}) to obtain an (the appropriate 'independent') subset Ia- Second, for 
a particular spacetime, we omit any scalars from Xa that can be obtained from syzygies 
defining that particular spacetime (e.g., defining the algebraic type of the spacetime, such 
as the Segre type or the Petrov type). For example, for a Ricci tensor corresponding to 
the algebraic form of a perfect fluid we could omit {R2,R3} (relative to {R,R1}). We 
consequently obtain the subset Xsa- e.g., 



Isa = {R, Rl, ■ ■ ■ , C 2 , . . .}. 

For the spacetimes under consideration the microgeometry is then completely characterized 
by the (sub)set of scalar curvature invariants Xsa- L 1 ! 

We now construct the new ordered set of functions Xsa by averaging the various scalar 
invariants of Xsa'- 

16 Including, for example, the differential Bianchi identities. 

17 This subset is not unique; perhaps a scheme should be adopted whereby higher order scalars should 
always be omitted; e.g., R 2 should be omitted relative to R. 
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Isa = {R, Rl, • • • , C 2 , . . .} 

where all of the original scalar invariants omitted from the original set X are replaced by a new 
set of functions obeying exactly the same algebraic properties (or syzygies) as Isa- Therefore, 
it is assumed that Isa comes equiped with these syzygies, so that we could construct the 
corresponding set I consisting of the members of Isa and all of the corresponding syzygies. 
Consequently, the set Isa is an ordered set of functions (scalar curvature invariants) on M 
which uniquely determines the macrogeometry with exactly the same algebraic properties as 
the original microgeometry. 

We shall refer to this proposal to obtain the set Isa and the associated averaged (macro) 
geometry as the Scalar Averaging Procedure. 

Note that in the cosmological application, in which the Ricci tensor is the fundemental 
object, the simple set Isa = {R, Rl} completely characterizes the averaged Ricci tensor (cf. 
[8]). However, this set does not completely determine the macrogeometry. In addition, we 
still need to know the trajectories along which the cosmological matter and null rays move; 
that is, macrogeometric effects not arising from the Ricci tensor alone. 

4.2 Cosmological model 

In the case of a cosmological model, we only need to be able to average the Ricci tensor (or 
Einstein tensor) that appears in the EFE. 
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4.2.1 Ricci tensor 



Since from C3 we have an effective set of field equations: 

R\ + C a b = T% (4.3) 

we only need to consider the macrogeometric Ricci tensor R a b here (the correlation tensor is 
obtained from the averaging procedure). The microgeometric Ricci tensor R a b is completely 
characterized by a set of scalar curvature invariants Xr. Averaging these scalar curvature 
invariants we obtain the set X^, which completely characterizes the macrogeometric Ricci 
tensor R a b . Since constructing the Ricci tensor from a set of scalar curvature invariants X^ is 
relatively simple compared to the corresponding problem for the Riemann tensor, and since 
the reduced set of scalar curvature invariants Xr is considerably smaller than X, we have 
considerably reduced the complexity of the problem in this new averaging approach. Indeed, 
for a Ricci tensor of the algebraic form of a perfect fluid, there are effectively (only) two 
independent scalar invariants, the Ricci scalar and a single Ricci eigenvalue (corresponding 
to the effective energy density, p, and pressure, p, of the perfect fluid). Therefore, in the 
context of the scalar averaging procedure, we have the set {R, Rl}. 

We also note that the syzygies of the macrogeometry Ricci tensor must be consistent with 
the syzygies of the energy-momentum tensor through the macroscopic field equation. If the 
energy-momentum tensor tensor is of the form of a perfect fluid; i.e., T a b = diag[—p,p,p,p], 
then T\ obeys a number of syzygies. A relation (equation of state) between p and p would 
be represented by a further algebraic syzygy. For example, defining A a b = T a b — ^TS a b , where 
T is the trace of T\ (compare with the form of the Ricci tensor in the EFE), the equation of 
state p + 3p = corresponds to the syzygy A 2 = ^A a b A b a , where A is the trace of A a b . o 

18 Since this particular syzygy can only be satisfied if p + 3p = or 5p + 3p = 0, if the cosmological matter 

satisfies appropriate energy conditions, this syzygy uniquely specifies the equation of state p + 3p = . 

19 Therefore, note that in order to obtain a macroscopic 'spatial curvature' it is necessary to include a 

cosmological pressure. 
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Finally, we note that if every member of Zr is zero (which uniquely characterizes a Ricci 
flat spacetime), then every member of X^, is also zero, and the macrogeometric Ricci tensor is 
also zero (flat). Therefore it follows that the average of a microvacuum R a b = gives rise to a 
macrovacuum R a b = (and the corresponding correlation tensor is zero, and the macroscopic 
field equations are trivial). Therefore, within this approach, no new effects due to averaging 
occur at the level of the macroscopic field equations when there is no cosmological (average) 
matter content (i.e., a non zero energy momentum tensor T a b is needed), o As we noted 
above, the primary focus of this paper is C3 and the effective macroscopic field equations. 

However, this does not imply that there are no averaging affects in the case of vacuum. 
Since from C4 we also need to be able to relate averaged quantities with observables and 
we need to know the trajectories along which the cosmological matter and null rays move, 
macrogeometric effects (not arising from the Ricci tensor alone) will indeed play a role. 
However, this is beyond the scope of the present work. It may be possible to express the 
effects in terms of kinematic scalars which are invariant with respect to all coordinate changes 
that preserve the fundamental congruence. 



4.2.2 Interpretation 



Finally, it is necessary to determine whether the correlations due to averaging alter the 
geometry or affect the effective energy-momentum tensor. This is partly a question of inter- 
pretation. This must be done within the context of the underlying cosmological model. We 
shall address this in the simple example in the next section. 

20 Note that in this approach, if the microgeometry is Ricci flat then the macrogeometry is Ricci flat. In 
vacuum, there will be no local variations (local inhomogeneities) in the Ricci tensor since there is local source 
with local inhomogeneities with a physical intrinsic scale. Indeed, the only inhomogeneities in the geometry 
must therefore come from (not from the energy-momentum tensor through the EFE but) physical boundary 
conditions, which would consequently be large scale inhomogeneities. Therefore, there would be no physical 
mechanism to introduce local inhomogeneities whose average would lead to a non-Ricci flat macrogeometry. 
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In particular, in the cosmological application it may be appropriate to reinterpret the 
averaging correlations as corrections to the matter fields through the EFE. Thus, concentrat- 
ing on the Ricci tensor again, writing G a b = T a b and defining Q\ =G a b — C a h , we can rewrite 
this as Q\ = T a b (by absorbing the correlation tensor C\ into T b ). We can now demand 
that the new macro-Ricci tensor 1Z a b (corresponding to Q a b ) has exactly the same algebraic 
properties as the Ricci tensor of the microgeomety, so that the averaging correlations are 
interpreted in the cosmological context as appropriate corrections to the matter fields. That 
is, we have defined the new averaged Ricci tensor TZ a b , derivable from some appropriate av- 
eraged geometry with an appropriate micro-metric, and some of the correction terms have 
been included in the new effective energy-momentum tensor T b . 
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5 Example: Static Spherically Symmetric Perfect Fluid 
Spacetimes 



We shall consider the specific example of a static spherically symmetric perfect fluid space- 
time in this section. This is an appropriate simple model for illustration since it can include 
an arbitrary function of one variable, there is a non- vanishing pressure, the averaging region 
does not change with time and there are no gravitational waves. 



Let us write the static spherically symmetric perfect fluid line element in the form 

ds 2 = -e f{r) dt 2 + e 9(r) [dr 2 + r 2 (d6 2 + sin 2 Odtf)]. (5.4) 



The two arbitrary functions f(r) and g{r) satisfy the differential constraint 



d 2 g ldg 1 df df dg if dg 



+ 



1 (df 



+ 



d 2 f 



dr 2 r dr r dr dr dr 2 \dr J 2 \dr ) dr 2 
in order for the Einstein tensor to be of the form of a perfect fluid. 



(5.5) 



The 16 CM polynomial scalar curvature invariants [30], as given in GRTensor [31], for 
the static spherically symmetric spacetime are then as follows (after the constraint has been 
applied). The 4 (linear, quadratic, cubic and quartic) Ricci tensor scalar invariants (which 
are related to the Ricci scalar and the three Ricci tensor eigenvalues) are: 



dr 2 \dr J \dr I \dr I \dr I r dr r dr 



(5.6) 



Rl 
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d 2 f 
dr 2 



df\ ( dg" 
dr I \dr 



dg" 
dr 



+ 



'df 
dr 



Adg 
r dr 



and R2 and R3 are proportional to (i?l) 3 / 2 and (-R1) 2 , respectively. 



(5.7) 



The 4 polynomial Weyl scalar invariants are: 



wm=-e 



-2g 



dr 2 



'df 
dr 



df\ ( dg" 
dr I \ dr 



2df 
r dr 



(5i 
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W2R oc (W1R) 3/2 and W1I = W2I = 0, and the mixed invariants are M3 = M2R oc 
(R1)(W1R), MA oc (R1)(W1R) 3/2 , W5R oc (R1) 3/2 (W1R) and MIR = Mil = M2I = 
M5I = 0. Note that GRTensor is also able to compute the differential scalar invariants 
diWeyl and <izRicci. 

In the case of a static spherically symmetric constant curvature spacetime with f(r) = 1, 
g(r) = (1 + kr 2 )~ 2 (k constant), we have that: 

R = 2Ak, Rl = 12k 2 (5.9) 

(R2 = 2Ak 3 , R3 = 8Ak 2 ), and all other scalar invariants (including the differential invariants) 
vanish. 

In the case of the static Schwartzchild-deSitter (Kottler [32]) spacetime with f(r) = 

1 - 2j f ~ = [g(r)}-\ we have that: 

(3?Ti^ — Qjtv^ 111 ^ 

R = 4A, W1R = —-, W2R = — —, diWeyl = -240(6m - 3r + Ar 3 )— (5.10) 

and all other scalar invariants (including diRicci) are zero. For the Schwarzchild vacuum 
solution A = 0. In the case of an Einstein space (with m — 0), the only non- vanishing scalar 
invariant is R = 4A. 

The approach outlined in the previous section is then as follows. Averaging the scalar in- 
variants we obtain, for example, (R) and (Rl) ((W1R) and so on, averaging eqns. fl5.6l) - fl5.8p . 
where (R) = R) together with an appropriately averaged constraint eqn. (|5.5|) . These are 
then the scalar curvature invariants of the averaged static spherically symmetric geometry. 
As noted above, we will focus on (R) and (Rl) here. 

5.1 Averaging Scales 

Suppose that the averaging scale is £ = L. Suppose also that there are inhomogeneities with 
scales [3~ l and A -1 , where of order unity (/3 ~ 1, L < i.e., global inhomogeneities) 
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and A" 1 is a small scale (A -1 << 1, A = where a is large (integer) so that X/L < 1; i.e., 
local small scale inhomogeneities). Hence, the function / varies slowly with respect to f3r, 
and varies quickly with respect to Ar (on smaller scales). 



Let us write the scale dependence explicitly as 

f = f((3r,\r), (5.11) 

where we effectively treat Ar as a separate variable. Consider r = r , and a neighbourhood 
of r , I = (r — -§,r + -§) of length L (the averaging region). Define r + r' £ I, so that 
r' £ (— f , f ) which parameterizes points in /. We can write (set (3=1) 

f(r) = f(r ,\(r + r')) (5.12) 

in J, where f(r ) = f(r , Ar ). Therefore, in J, f(r) has small scale variations with respect 
to Ar'. We can write 

f(r)~f(r ,\r )+O(n) (5.13) 
where \i is the small scale of the amplitude of inhomogeneities in I (/i « 1). 

We now average f(r) over these small scale inhomogeneities in /, and define the average 
(f(r )) = (f(r,\r))\ r=ro by: 

1 r L / 2 

(f(ro)) = j f(r ,X(r Q + r')dr' (5.14) 

L J-L/2 

(i.e., we are effectively 'averaging over Ar'). Notice that there are a number of scales in the 
problem: j « 1, pi « 1, ^ = j < 1. Note also that 

where the second term is a dominant fast varying term. 

We now Fourier analyse the functions f(r) and g(r) in the averaging region / with respect 
to (Ar') [A- 1 = a/L}: 

fir) = f(r ) + fif^ro) sin {^f^j + E ^fn(r ) sin (^JL^j 
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((X7TT \ / TlOtlTT \ 

~L~\ + P n 9n{r Q ) sin I ^ j (5.16) 

where /i(r ), f n (fo), 9l{ r o), <7n(?"o) are slowly varying functions of r (i.e., their derivatives 
are small). Note that 

</(ro)> = /(r ) 

(<?(r )> = ^(r ). (5.17) 



Calculating, we obtain (with an abuse of notation) 

au ( aitr'\ fi 2 a ( 2a7ir'\ 

9r{r) = g r {ro) + -^^91 cos I — — I + —2icg 2 cos I — — I 

. f airr'\ 9 , ( 2a-nr'\ _ f u 2 a 2 \ 

+lig[ sin _ + fi 2 g' 2 sin f — — J + 0{^) + O ^ 

/ia 2 / a7rr'\ u 2 a 2 / 2anr'\ 

9rr{r) = g r r{r ) - ~JY n 9i sm I — — I - -jyAng 2 sin I — — I 

+0(^) + 0(^j (5.18) 

We can now average the constraint equation (I5.5P [using various trigonometric formulae in 
evaluating the integrals]: To zeroth order we obtain the constraint in terms of f(ro) and 
g(ro) [effectively eqn. (15.51) ]. To next leading order we obtain 

h9i + \g\ - \fl = (5.19) 

(and to higher orders: f[g[ + \(g[) 2 - ±(f[) 2 = 0, f 2 g 2 + \{gl) - \f% = 0, etc.). For example, 
using the first of eqns. (15 . 1 8[) to expand g r (r) and f r (j) and integrating we obtain 

I df dg\ . 7r 2 (a\ 2 2 1 , , 



dr dr 



(5.20) 



fr(r )gr(r ) + y/ifi-i (^j ^ + ^fWi^ 2 

+ Y h92 (z) 2/i4 + l fWy + i^f) + ^ 6) ' 

Therefore, by a direct calculation (and using the constraint (15.191) . we find that 

(R) = R(ro) (l + - ^e-(-) (j)\ 2 \ (fin + O (I)) + (5.21) 
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We also find that 

(Rl) = RHro) |l + O (£) } + ^fF(r ) (| ) V~(/i0l), (5.22) 
where F(r Q ) is a specific function of r (not explicitly displayed here). 

5.2 Interpretation 

The O(t-) correction in the first term of (R) in eqn. (I5.2ip is a higher order renormalization 
term. To lowest order we have that 

(R) = R(ro) - ^e-(-) {hgi) 

and hence 

= i? 2 (r ) - ^G(r )e (5.23) 

where G(r ) = e~ 9 ^°)i?(r ) and e = (/i<7i) — constant. Also, we have that 

(Rl) = Rl(r ) + y#(r )e, (5.24) 

where H(r ) = Rl(r )F(r ). For a space of constant curvature k with R(r ) = 24k, Rl(r ) = 
12k 2 = ^(24/c) 2 , we have that G(r ) = 7687 fc 2 , H(r Q ) = — ^^ 2 7o (where 7 is a constant). 
Therefore, we have 

(R) 2 = R 2 (r ) + 2(24A;) 2 (-7oe) 

(J21) = m(r ) + l(24/c) 2 (- 7 oe). (5.25) 

and we can interpret the average correlations as contributing a small, 0(a 2 n 2 L~ 3 ), constant 
curvature term, arising from the averaging of local inhomogeneities in the micro-Ricci tensor, 
to the smooth macro- Ricci tensor (consistent with the results of [14J). 
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